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The properties of multiple stellar systems have long provided important empirical con- 
straints for star formation theories, enabling (along with several other lines of evidence) a 
concrete, qualitative picture of the birth and early evolution of normal stars. At very low 
masses (VLM; M < 0. 1 Mq), down to and below the hydrogen burning minimum mass, our 
\ understanding of formation processes is not as clear, with several competing theories now under 

consideration. One means of testing these theories is through the empirical characterization 
of VLM multiple systems. Here, we review the results of various VLM multiplicity studies to 
date. These systems can be generally characterized as closely separated (93% have projected 
^ ■ separations A < 20 AU) and near equal-mass (77% have M2/M1 > 0.8) occurring infrequently 

^ ' (perhaps 10-30%). Both the frequency and maximum separation of stellar and brown dwarf 

, binaries steadily decrease for lower system masses, suggesting that VLM binary formation 

' and/or evolution may be a mass-dependent process. There is evidence for a fairly rapid decline 

^ I in the number of loosely-bound systems below ~0.3 Mq, corresponding to a factor of 10-20 

• 1^ . increase in the minimum binding energy of VLM binaries as compared to more massive 

' stellar binaries. This wide-separation "desert" is present among both field (~l-5 Gyr) and 

JH , older (> 100 Myr) cluster systems, while the youngest (<10 Myr) VLM binaries, particularly 

- - 1 those in nearby, low-density star forming regions, appear to have somewhat different systemic 

properties. We compare these empirical trends to predictions laid out by current formation 
theories, and outline future observational studies needed to probe the full parameter space of 
the lowest mass multiple systems. 
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1. INTRODUCTION 

The frequency of multiple systems and their properties 
are key constraints for studies of stellar formation and evo- 
lution. Binary and multiple stars are common in the Galaxy, 
and the physical properties of the components in these sys- 
tems can be significantly influenced by dynamical and co- 
evolutionary processes. Furthermore, successful theories of 
star formation must take into account the creation of multi- 
ples and empirical multiplicity trends as functions of mass, 
age and metallicity. 



The main focus of this review is multiplicity in very low 
mass (VLM; M < 0.1 M0) stars and brown dwarfs. How- 
ever, to put these results in the proper context, we start with 
a brief review of our current understanding of multiplic- 
ity among higher mass stars (also see chapter by Duchene 
et al.). The standard references for binary frequency are 
Duquennoy & Mayor (1991, hereafter DM91; also Abt and 
Levy, 1976; Abt, 1978; Mayor et al, 1992) for solar-type 
stars and Fischer and Marcy (1992, hereafter FM92; also 
Henry and McCarthy, 1990; Reid and Gizis, 1997a; Halb- 
wachs et al., 2003; Delfosse et al., 2004) for early-type M 
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dwarfs. The DM91 survey combined spectroscopic, astro- 
metric and direct imaging of 164 G dwarfs; 44% of those 
stars were identified as binaries, with incompleteness cor- 
rections increasing the binary fraction to ftin ~ 65%. 
These corrections include 8% attributed to VLM compan- 
ions; as discussed further below, more recent observations 
show that the actual correction is much lower. The FM92 
survey covered 72 M2-M5 dwarfs within 20 parsecs, and 
derived fbm = 42 ± 9%, significantly lower than the DM9 1 
G-dwarf survey. While both surveys include nearby stars, 
neither comprises a volume-complete sample. 
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Fig. 1. — Mass and mass ratio distributions of companions to 
late-F to K-type dwarfs within 25 parsecs of the Sun {Reid et al., 
2002a), segregated by projected separation/orbital semimajor axis. 
The left panels plot the mass distribution of companions, with the 
dashed lines providing a schematic representation of the initial 
mass function. The right panels plot the mass-ratio distributions 
(dotted histograms), with the solid histogram showing the mass 
ratio distribution for VLM dwarfs (no segregation of separations; 
see Fig. 3). These distributions are normaUzed at the g = 1 bin. 

Recent surveys of solar-type stars have concentrated on 
VLM companions. Radial velocity (RV) surveys (e.g., 
Marcy and Butler, 2000; Udry et al, 2003) have shown 
that less than 0.5% of solar-type stars have brown dwarf 
companions within ^ 5 AU. Guenther et al. (2005) find 
Ibin < 2^ projected separations A < 8 AU among 
Hyades stars; this is in contrast with f^in ~ 13% for stellar- 
mass companions at those separations (DM91). At larger 
separations, imaging surveys of young Solar neighborhood 
stars (members of the TW Hydrae, Tucanae, Horologium 
and /3 Pic associations; Zuckerman and Song, 2004) find 
fun ~ 6 ± 4% for A > 50 AU (Neuhauser et al. 
2003), similar to the brown dwarf companion fraction 
measured for field stars for separations of 30-1600 AU 
(Metchev, 2005). These fractions are ~3 times lower than 



the hydrogen-burning companion rate over the same sepa- 
ration range. At the widest separations (A > 1000 AU), 
Gizis et al. (2001) find that solar- type stars have compa- 
rable numbers of brown dwarf and M dwarf companions, 
although this result is based on a very small number of 
VLM companions. 

Besides the overall binary fraction, the mass distribution 
of companions sets constraints on formation models. Fig. 1 
shows the results for late-F to K stars (0.5 < {B — V) < 
1.0) within 25 parsecs of the Sun, breaking down the sample 
by projected separation/orbital semi-major axis. The left 
panels compare the mass distribution of companions against 
a schematic representation of the initial mass function (Reid 
et al, 1999, 2002a); the right panels compare the mass ra- 
tio {q = M2/M1) distributions against the VLM dwarf data 
assembled in this review (cf.. Fig. 3). Clearly, low q bi- 
nary systems are more common at all separations among 
solar- type stars than in VLM dwarfs. We return to this is- 
sue in Section 2.2.3. At small separations (A < 10 AU), 
there is an obvious deficit of low-mass companions (with 
the exception of planetary companions) as compared to the 
distribution expected for random selection from the field- 
star mass function. The notorious brown dwarf desert (e.g., 
Marcy and Butler, 2000) extends well into the M dwarf 
regime. This result is consistent with the original analysis 
of Mazeh and Goldberg (1992) of the mass ratio distribution 
of spectroscopic binaries, although their more recent study 
of proper motion stars {Goldberg et al, 2003) finds a bi- 
modal distribution, with peaks at q ^0.8 and ^0.2 (see also 
Halbwachs et al., 2003). The deficit in low-mass compan- 
ions is less pronounced at intermediate separations, while it 
is possible that observational selection effects (e.g., sensi- 
tivity limitations) might account for the small discrepancy 
for g < 0.2 in the wide-binary sample. 

In the case of M dwarfs, attention has focused on the 
nearest stars. Delfosse et al. (2004) recently completed a 
spectroscopic and adaptive optics (AO) imaging survey of 
M dwarfs within 9 parsecs that is effectively complete for 
stellar mass companions. Combining their results with the 
imaging surveys by Oppenheimer et al. (2001) and Hinz et 
al. (2002), they derive an overall binary fraction of 26% 
for M dwarfs. For a more detailed breakdown with spec- 
tral type, we can turn to the northern 8-parsec sample {Reid 
and Gizis, 1997a; Reid et al, 2003). Those data indicate 
binary fractions of 24+^^% for spectral types M0-M2.5 
(4/17 systems), 2lt\% for M3-M4.5 (12/45 systems) and 
3llg'^% for M5-M9 (5/16 systems; uncertainties assume a 
binomial distribution), where the spectral type refers to the 
primary star in the system; the overall binary frequency is 
fbin = 27^4%. These results, based on volume-limited 
samples, confirm that M dwarfs have significantly lower 
multiplicity than more massive solar-type stars. ^ This is 



Even with 30% binarity for M dwarfs, most stars still reside in multiple 
systems. As a numerical example, consider a volume-limited sample of 
100 stellar systems: 20 are type G or earlier, 10 are type K and 70 are type 
M. Assuming binary fractions of 70%, 50% and 30%, respectively, these 



2 



consistent with an overall trend of decreasing multiplicity 
with decreasing mass (cf., A- and B-stars have overall mul- 
tiphcity fractions as high as 80%; Shatsky and Tokovinin, 
2002; Kouwenhoven et ah, 2005). These changes in multi- 
plicity properties with mass among hydrogen-burning stars 
emphasize that we must consider VLM dwarfs as part of a 
continuum, not as a distinct species unto themselves. 

2. OBSERVATIONS OF VERY LOW MASS BINA- 
RIES 

2.1 Very Low Mass Binary Systems 

With the discovery of hundreds of VLM dwarf stars and 
brown dwarfs over the past decade (see reviews by Basri, 
2000; Oppenheimer et ah, 2000; and Kirkpatrick, 2005), it 
is now possible to examine systems with primaries down to 
100 times less massive than the Sun. In this regime, forma- 
tion mechanisms are under considerable debate (see chap- 
ters by Bonnell et al., Goodwin et al., Klein et al., Luhman 
et al., and Whitworth et al.). Hence, accurate assessment of 
the multiplicity and systemic properties of VLM stars and 
brown dwarfs are essential for constraining current theoret- 
ical work. 

Searches for VLM binaries — defined here as having a 
total system mass M-tot < 0.2 M© and primary mass Mi 
< 0.1 Mq (cf., Siegler et al, 2005) — have been con- 
ducted predominantly through high resolution imaging sur- 
veys, using both ground-based (including natural and, quite 
recently, laser guide star adaptive optics [AO]) and space- 
based facilities. Major surveys have targetted both nearby 
field sources {Koerner et al, 1999; Reid et al., 2001; Bouy 
etal, 2003; Burgasser et al, 2003; Close etal, 2002, 2003; 
Gizis et al, 2003; Siegler et al, 2003, 2005; Law et al., 
2006; Allen et al. in preparation; Billeres et al. in prepara- 
tion; Burgasser et al. in preparation; Reid et al in prepa- 
ration) and young clusters and associations {Martin et al, 
1998, 2000a, 2003; NeuhaUser et al, 2002; Kraus et al, 
2005; Luhman et al, 2005; Bouy et al, 2006). A smaller 
number of high resolution spectroscopic surveys for closely 
separated binaries have also taken place {Basri and Martin, 
1999; Joergens and Guenther, 2001; Reid et al, 2002b; 
Guenther and Wuchterl, 2003; Kenyan et al, 2005; Joer- 
gens, 2006). Only one echpsing system has been discovered 
so far via photometric monitoring {Stassun et al, 2006). 
Observations leading to the identification of low mass mul- 
tiple systems has been accompanied by resolved photome- 
try and spectroscopy, allowing characterization of the col- 
ors, luminosities and spectral characteristics of several bi- 
nary components. Astrometric and radial velocity monitor- 
ing has lead to mass measurements or constraints for five 
VLM systems to date {Basri and Martin, 1999; Lane et al, 
2001b; Bouy et al, 2004a, Brandner et al, 2004; Zapatero 
Osorio et al . 2004: Stassun et al, 2006). 



100 systems include 140 stars, 80 in binaries and 60 in isolated systems. 
Higher order multiples only serves to increase the companion fraction. 



In Table 1 we list 75 VLM binary systems published 
in the Uterature or reported to us as of 2005. The mass 
criteria correspond to field dwarf binary components later 
than spectral type '~M6; younger systems may include ear- 
lier spectral types. Table 1 provides a subset of the com- 
piled data for these sources, given in more complete detail 
through an onhne database maintained by N. Siegler (See 
http://papercUp.as.arizona.edu/~nsieglerA^LMJbinaries.). 

2.2 General Properties of VLM Binaries 

Large-scale, high resolution imaging surveys in the field 
have converged to similar conclusions on the general prop- 
erties of VLM field binaries. Compared to their higher mass 
stellar counterparts, VLM binaries are 

• rarer {fbin « 10 — 30%; however, see discussion be- 
low); 

• more closely separated (93% have A < 20 AU); 

• and more frequently in near-equal mass configura- 
tions (77% have q > 0.8). 

Analogous imaging surveys in young open clusters (e.g., 
Pleiades, a Persei) find similar trends, although the youngest 
(<10 Myr) associations (e.g., Chamaeleon I, Upper Scor- 
pius, Orion) appear to exhibit somewhat different prop- 
erties. We discuss these broad characterizations in detail 
below. 

2.2.1 The Binary Fraction 

Magnitude-hmited imaging surveys for VLM stars and 
brown dwarfs in the field with spectral types M6 and later 
have generally yielded observed binary fractions of ^20%; 
taking into consideration selection effects (e.g., Burgasser 
efflZ., 2003) lowers this fraction to 7-15% for A > 2-3AU 
and g > 0.4-0.5 {Bouy et al, 2003; Burgasser et al 2003; 
Close et al, 2003; Siegler et al, 2005). Burgasser et 
al (2003) deduced fbin = for a small sample of 

L and T dwarfs using the 1 /Vmax technique {Schmidt et al , 
1968); Bouy et al (2003) deduced a volume-limited frac- 
tion of fbin ^ 15%. Over the same separation (A > 2 AU) 
and mass ratio {q > 0.5) ranges, these multiplicity rates are 
less than half of those of M dwarfs (FM92; Close et al, 
2003) and G dwarfs (DM91; Bouy et al, 2003). Similarly, 
HST imaging surveys of the 125 Myr Pleiades open clus- 
ter {Martin et al, 2000a, 2003; Bouy et al, 2006) found a 
resolved binary fraction of 13-15% for A > 7 AU for com- 
ponents at and below the hydrogen burning limit. On the 
other hand, Kraus et al (2005) found fbin = 25^1^% for a 
small sample of 0.04-0. 1 Mq members of Upper Scorpius 
over the range A = 5-18 AU, somewhat higher than, but 
still consistent with, other field and open cluster results. 

One problem with resolved imaging surveys is their 
inherent selection against tightly bound systems (A < 
2 — 3 AU for the field dwarfs and nearby associations, 
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A < 10 — 15 AU for more distant star forming regions). 
Here, one must generally turn to high resolution spectro- 
scopic surveys of VLM stars, currently few in number and 
with as yet limited follow-up. Reid et al. (2002b) de- 
duced a double-lined spectroscopic binary (SB2) fraction of 
6^2% for a sample of M7-M9.5 field dwarfs. Guenther and 
Wuchterl (2003) identified two SB2s and marginally sig- 
nificant RV variations in the active M9 LP 944-20 (which 
they attribute to either the presence of a low-mass compan- 
ion or magnetic-induced activity) in a sample of 25 M5.5- 
L1.5 field and cluster dwarfs. Including all three objects 
implies an observed binary fraction of 12^4°%, although 
this value does not take into consideration selection biases. 
Joergens (2006) detected one RV variable, the M6.5 Cha 
Ha8, among a sample of 9 VLM stars and brown dwarfs in 
the 2 Myr Cha I association, implying an observed fraction 
of 1 1^4*%, again subject to sampling and selection biases. 
Kenyan et al. (2005) identified four possible spectroscopic 
binaries (SBs) among VLM stars and brown dwarfs in the 
3-7 Myr a Orionis cluster on the basis of RV variations 
over two nights. They derive > 7 — 17% for A < 
1 AU (after correcting for selection effects) and a best-fit 
fraction of 7-19% (for their Sample A) depending on the 
assumed underlying separation distribution. However, none 
of the sources from this particular study have had sufficient 
follow up to verify RV variability, and cluster membership 
for some of the targets have been called into question. A 
more thorough analysis of sensitivity and sampling biases 
in these SB studies has been done by Maxted and Jeffries 
(2005), who find fun = 17-30% for A < 2.6 AU, and an 
overall binary fraction of 32^5% (assuming fbin = 15% 
for A > 2.6 AU). This result suggests that imaging stud- 
ies may be missing a significant fraction of VLM systems 
hiding in tightly- separated pairs. However, as orbital prop- 
erties have only been determined for two SB systems so far 
(PPl 15, Basri and Martin, 1999; and 2MASS 0535-0546, 
Stassun et al., 2006), individual separations and mass ratios 
for most VLM SB binaries remain largely unconstrained. 

Two recent studies {P infield et al., 2003; Chapelle et al, 
2005) have examined the fraction of unresolved (overlu- 
minous) binary candidates among VLM stars and brown 
dwarfs in young associations. Contrary to other studies, 
these groups find much larger binary fractions, as high as 
50% in the Pinfield et al. study of the Pleiades and Prae- 
sepe for q > 0.65. This study also finds a binary fraction 
that increases with decreasing mass, in disagreement with 
results in the field (see below); the Chapelle et al. study 
finds evidence for the opposite effect in the 0.9 Gyr Prae- 
sepe cluster. Both studies have been controversial due to 
the lack of membership confirmation, and hence likelihood 
of contamination; and the possible influence of variability 
on the identification of overluminous sources. Neverthe- 
less, both of these studies and the SB results suggest that a 
higher VLM binary fraction than that inferred from imaging 
studies, perhaps 30% or more, is possible. 

2.2.2 The Separation Distribution 



Fig. 2 plots the histogram of projected separations/orbital 
semimajor axes for 70 binaries in Table 1 (SB systems with- 
out orbital measurements are not included). This distribu- 
tion exhibits a clear peak around 3-10 AU, with 53±6% 
of known VLM binaries encompassing this range. Again, 
because imaging surveys (from which most of the objects 
in Table 1 are drawn) can only resolve systems down to a 
minimum angular scale (typically 0'.'05 — O'.'l for AO and 
HST programs), the decUne in this distribution at small 
separations is Ukely a selection effect. Results from SB 
studies remain as yet unclear in this regime. Basri and 
Martin (1999) have suggested that very close binaries are 
common based on the detection of one (PPl 15) in a small 
spectral sample The analysis of Maxted and Jeffries (2005) 
suggest that there may be as many or more binaries with 
A < 3 AU as fliose with A > 3 AU. At the extreme. Pin- 
field et al. (2003) estimate that 70-80% of VLM binaries in 
the Pleiades have A < 1 AU, although this result has not 
been corroborated by similar studies in the Pleiades {Bouy 
et al, 2006) and Praesepe {Chapelle et ah, 2005). In any 
case, as the peak of the observed separation distribution 
lies adjacent to the incompleteness Umit, closely separated 
systems Ukely comprise a non-negUgible fraction of VLM 
binaries. 
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Fig. 2. — Distribution of separations/orbital semimajor axes 
for known VLM binary systems (Table 1). The number 
of VLM binary systems in each 0.3 dex bin is labelled, 
and uncertainties (vertical lines) are derived from a bino- 
mial distribution. Note that SBs with unknown separations 
are not plotted but included in the total number of bina- 
ries for scaling the distribution. The distribution peaks at 
A ~ 3 — 10 AU, with steep dec fines at shorter and longer 
separations. While there is Ukely observational incomplete- 
ness for A < 3 AU, the sharp drop in binary systems with 
A > 20 AU is a real, statistically robust feature. The shaded 
bins represent the 8 systems with ages < 10 Myr. While the 
statistics are still small, the separation distribution of these 
young binaries is flatter, and suggests a peak at wider sepa- 
rations than that of the field and older cluster binaries. 

The steep decline in the separation distribution at larger 
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separations is, on the other hand, a statistically robust fea- 
ture. While high resolution imaging surveys are limited in 
this domain by field of view (typically 10-20" for HST and 
AO studies), this only excludes systems with A > 150 AU 
for a typical VLM field source (distances ~ 30 pc) or 
A > 200 — 1000 AU for young cluster systems. Even 
wider separations for hundreds of VLM field dwarfs should 
be detectable - and are not found - in the original sur- 
veys from which they were identified (e.g., 2MASS, DE- 
NIS and SDSS; however, see Billdres et al, 2005). In open 
clusters, deep imaging has demonstrated a consistent lack 
of wide companions to VLM dwarfs. An upper hmit of 
fun < 8% for A > 11 AU is derived for the 90 Myr a Per 
open cluster {Martin et ah, 2003), similar to the 5% upper 
hmit for A > 15 AU measured for 32 VLM members of 
the 2 Myr IC 348 cluster {Luhman et al, 2005). Lucas et 
al. (2005) measure an upper limit of 2% for wide VLM bi- 
naries (A > 150 AU) in the 1 Myr Trapezium cluster based 
on a two-point correlation function. In contrast, 93% of the 
known VLM binaries have A < 20 AU. Hence, a "wide 
brown dwarf binary desert" is evidenced for VLM stars and 
brown dwarfs {Martin et al. , 2000a), a potential clue to their 
formation. 

While survey results have generally been negative 
for wide VLM binaries, two — 2MASS J11011926- 
7732383AB {Luhman, 2004; hereafter 2MASS 1101- 
7732AB) and DENIS J055146.0-443412.2AB {Billeres et 
al, 2005, hereafter DENIS 0551-4434AB) - have been 
identified serendipitously. These systems have projected 
separations > 200 AU, over 10 times wider than the vast 
majority of VLM binary systems. A third low mass binary 
not included in Table 1, GG Tau BaBb (a.k.a. GG Tau/c; 
Leinert et al, 1991; White et al, 1999), with estimated 
primary and total system masses of 0.12 and 0.16 M©, 
respectively, also has a projected separation greater than 
200 AU. Interestingly, two of these three systems are mem- 
bers of very young, loose associations. We discuss these 
source further in §2.4.2. 

The separation distribution of VLM stars therefore peaks 
at or below ~3-10 AU, corresponding to orbital periods of 
<40 yr. This is quite different from the separation distribu- 
tion of G dwarfs, which shows a broad peak around 30 AU 
(periods of ~170 yr; DM91); and the M dwarf distribu- 
tion, which peaks between 4-30 AU (periods of 9-270 yr; 
FM92). There is a suggestion in this trend of decreasing 
separations as a function of mass, as discussed further be- 
low. 

2.2.3 The Mass Ratio Distribution 

Fig. 3 shows the distribution of mass ratios for 70 of the 
binaries in Table 1 (not including SBs without mass esti- 
mates). These ratios were derived by a variety of methods, 
including comparison of component fluxes to evolutionary 
models (e.g., Chabrieretal, 2000), analytic relations (e.g.. 
Burrows et al, 2001) and direct estimates from orbital mo- 
tion measurements. Despite these different techniques, a 



comparison of all the data shows congruence with individ- 
ual studies. The mass ratio distribution for VLM systems is 
strongly peaked at near-unity ratios; over half of the known 
VLM binaries have q > 0.9 and 77^^% have q > 0.8. 

As with the separation distribution, it is important to con- 
sider selection effects in the observed mass ratios. Most per- 
tinent is the detectabiUty of secondaries in low q binaries, 
which may be too faint for direct imaging or of insufficient 
mass to induce a measureable RV variation in the primary's 
spectrum. The former case is an important issue for field 
binaries, as low mass substellar companions fade to obscu- 
rity over time. However, most imaging and spectroscopic 
surveys to date are sensitive down to g > 0.5, while a sharp 
dropoff is clearly evident at the highest mass ratios. Hence, 
while the number of low mass ratio systems may be under- 
estimated, the g ~ 1 peak is not the result of this bias. 
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Fig. 3. — Mass ratio distribution of known VLM binary 
systems (Table 1). The number of VLM binary systems 
in each 0.1 fractional bin is labelled, and uncertainties are 
derived from a binomial distribution. Note that SBs with 
unknown mass ratios are not plotted and not included in 
the total number of binaries when scaling the distribution. 
The distribution peaks near unity for binary systems with 
A > 3 — 4 AU, and matches a power law. Note that incom- 
pleteness is hkely for q < 0.6. The shaded bins represent 
the eight systems with ages < 10 Myr. While the statistics 
are still small, the mass ratio distribution of these young 
systems suggests a flatter distribution than that of field and 
older cluster binaries. 

A second effect is the preferential discovery of unre- 
solved equal-mass systems in wide field surveys. As such 
systems are twice as bright as their single counterparts, they 
are ^3 times more likely to be found than single sources in 
a magnitude-hmited survey. Systems with lower mass ra- 
tios are not as overluminous and less affected by this bias. 
Burgasser et al (2003) examined this impact of this bias on 
a small sample of L and T dwarf binaries and found it to 
be significant only for q < 0.6. Hence, this bias cannot be 
responsible for the q ^ 1 peak. 

VLM (field and open cluster) binaries therefore show a 
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clear preference for equal mass systems, in contrast to the 
majority of F-K stellar systems (Fig. 1). It is worth noting 
that M dwarfs in the 8 pc sample show a similar, although 
less pronounced, g ~ 1 peak {Reid and Gizis, 1997a), again 
suggesting a mass-dependent trend. 

2.2.4 Higher Order Multiples 

Thus far we have focused on VLM binaries, but higher 
order multiples (triples, quadruples, etc.) are also abun- 
dant among more massive stars, comprising perhaps 15- 
25% of all multiple stellar systems (Tokovinin, 2004; see 
chapter by Duchene et al.). Several VLM binaries are com- 
ponents of higher order multiple systems with more mas- 
sive stars. Burgasser et al. (2005a) have even suggested 
a higher binary fraction for brown dwarfs that are widely- 
separated companions to massive stars. Higher order mul- 
tiples are currently rare among purely VLM systems, how- 
ever. The LP 213-67/LP 213-68AB system is one excep- 
tion, with the three components (spectral types M6.5, MS 
and LO) forming a wide hierarchical triple with separations 
of 340 AU and 2.8 AU {Gizis et al, 2000a; Close et al, 
2003). DENIS 0205-1159AB may also have a third com- 
ponent, marginally resolved through high resolution imag- 
ing {Bouy et al, 2005). Considering both systems as VLM 
triples, the ratio of high-order multiples to binaries is only 
3^1%, quite low in comparison to higher mass stars. This 
may be due to selection effects, however, as the already 
tight separations of VLM binaries implies that the third 
component of a (stable) hierarchical triple must be squeezed 
into an extremely small orbit. Indeed, this could argue 
against a large fraction of higher order VLM systems. On 
the other hand, undiscovered wide tertiaries (as in LP 213- 
67/LP 213-68AB) may be present around some of these sys- 
tems. Additional observational work is needed to determine 
whether higher order VLM multiples are truly less common 
than their stellar counterparts. 

2.3 Statistical Analysis: Bayesian Modeling 

To examine the observed binary properties of resolved 
VLM stars in more detail, we performed a Bayesian sta- 
tistical analysis of imaging surveys to date. The Bayesian 
approach allows the incorporation of many disparate data 
sets, and the easy assimilation of non-detections, into a uni- 
fied analysis of a single problem (Sivia, 1996). We focused 
our analysis on the surveys of Koerner et al. (1999); Reid et 
al (2001); Bouy et al (2003); Close et al (2003); Gizis et 
al. (2003); Siegler et al. (2005); and Allen etal. (in prepara- 
tion). The Bayesian statistical method employed is similar 
to that described in Allen et al. (2005). 

We first constructed a set of parameterized companion 
distribution models in terms of orbital semi-major axis (a) 
and companion mass ratio. For the semi-major axis distri- 
bution we use a Gaussian in log AU given by: 

P(ao, aa) = ^_e-(i°s(«)-i°g(«o))V2a2 (1) 



where ao is the peak of the Gaussian and (Ja is the logarith- 
mic half-width, both variable parameters. This formulation 
is prompted by the results of DM91 and FM92 (however, 
see Maxted and Jeffries, 2005). For the mass-ratio model, 
we assume a power law of the form: 



where the normalization factor N is defined to be the overall 
binary fraction (i.e., /bin), and 7 is a variable parameter. 

In order to compare the model distributions to the data, 
we transform them to observables, namely the log of the 
projected separation (log A) and the difference in magni- 
tude between the secondary and the primary (AM). The 
former is computed by transforming the semi-major axis 
distribution as: 

A = a\^cos'^{(j))sin'^{i) + sin'^{(p), (3) 

where we assume a uniform distribution of circular orbits 
over all possible inchnations (i) and phases ((/>). The trans- 
formation of the q distribution to a AM distribution is done 
by assigning each mass ratio a range of possible luminosi- 
ties for ages between 10 Myr and 10 Gyr using evolutionary 
models from Burrows et al (2001). 

The transformed model distributions are then compared 
to the observed distributions via a Bayesian statistical 
method, as described in Allen et al (2005). The models are 
directly compared to the data after being convolved with a 
window function, which describes how many times a bin in 
observational space (A, AM) was observed in a particular 
survey. In this way we do not analyze the models where 
there is no data, and the relative frequency of observations 
is taken into account. 

The output posterior distribution is four dimensional 
(log(ao), Ua, N, 7) and is impossible to display in its en- 
tirety. Instead, we show marginahzed distributions (Fig. 4), 
collapsing the posterior distribution along different param- 
eter axes. These distributions have a non-neghgible dis- 
persion, as parameters spaces outside the observational 
window function (e.g., very tight binaries) add consider- 
able uncertainty to the statistical model. Nevertheless, 
the distributions are well-behaved and enable us to derive 
best-fit values and uncertainties for the various parame- 
ters. The overall binary fraction is reasonably well con- 
strained, N = 22^4%, with a long tail in its probability 
distribution to higher rates. The remaining parameters are 
log(ao) = OMt°oil ^og{AU), = Q.2A+J,°^^ \og{AU), 
and 7 = 4.8^ i g (all Usted uncertainties are 68% confidence 
level). 

The mass ratio and projected separation distributions in- 
ferred from the best- fit parameters are shown in Fig. 5. 
The best-fit binary fraction is 22%, but after applying our 
window function the expected resolved fraction is ~17%, 
sUghtly higher than but consistent with the observed fhin 
from imaging surveys (§ 2.2.1). The best-fit mass ratio dis- 
tribution is highly peaked near g = 1, similar to the data 
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Fig. 4. — Posterior probability distributions of the four 
companion model parameters: a) overall binary fraction 
(N); b) center of the semi-major axis distribution (log(ao)); 
c) width of the semi-major axis distribution ((To); d) mass 
ratio distribution power law index (7). 



but somewhat flatter than observed due to selection effects 

in the empirical samples. This nevertheless confirms that 
the mass ratio distribution is fundamentally peaked towards 
high q values. 

The best fit value for the peak of the semimajor axis dis- 
tribution is ~7 AU, implying a peak in the projected sepa- 
ration distribution of about 3.5 AU, matching well with the 
data (Fig. 5b). The best fit width of this distribution is quite 
narrow, implying very few wide systems (> 20 AU ~ 1%) 
and very few close systems (< 1 AU ~ 2-3%). It is impor- 
tant to stress that the imaging data provide weak constraints 
on closely-separated binaries, and the latter fraction may 
be somewhat higher (cf., Maxted and Jeffries, 2005). On 
the other hand, the constraint on the wide binary fraction 
(1% or less) is the most robust result of this analysis. Be- 
tween all of the surveys considered here there are over 250 
unique fields that have been probed for companions out to 
hundreds of AU with no detections. Hence, such pairings 
are exceptionally rare. 

2.4 Discussion 



2.4.1 On the Preference of Tight Binaries 

The sharp decline in the VLM binary fraction for A > 
20 AU is not a feature shared with more massive stellar sys- 
tems, which can extend from 0.1 AU to 0.1 pc. However, 
the decline is consistent with the observed trend of smaller 
mean separations, and smaller maximum separations, from 
A to M field binaries. This is demonstrated in Fig. 6, which 
plots projected separations/semimajoraxes versus total sys- 
tem mass for stellar and substellar field and cluster binaries. 
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Fig. 5. — (Top) The fraction of VLM binaries with a given 
q for the best fit model (solid line), the best fit model view 
through the window function (dashed line), and the data 
used in the Bayesian analysis (dot-dashed Une). Note how 
the window function over-emphasizes the high mass ratio 
systems. (Bottom) The projected separation distribution for 
the best fit model (lines are the same as the top panel). 

The maximum separations (Amax) of these systems show 
a striking dependence on total system mass. Prior to the 
discovery of the wide pairs 2MASS 1 101-7732AB and DE- 
NIS 0551-4434AB, Burgasser et al (2003) found a power- 
law relation between Amax and total system mass, Amax = 
14OO(M(ot/M0)2 AU, that appeared to fit all VLM systems 
known at that time. Similarly, Close et al. (2003) found 
a linear relation of A„„,,,. = 23.2(M,o,/0.185 Mg) AU for 
VLM binaries, corresponding to a minimum escape veloc- 
ity Vesc = 3.8 km s~^. This was greater than a minimum 
value of Vesc = 0.6 km s^^ inferred for more massive stellar 
systems, and both results indicate that lower mass binaries 
are progressively more tightly bound. Close et al. (2003) 
further pointed out a possible "break" in the minimum bind- 
ing energies of stellar and VLM binaries, also shown in 
Fig. 6. Around Miot ~ 0.3 M©, the majority of wide VLM 
systems appear to be 10-20 times more strongly bound than 
the widest stellar systems. 

More recently, exceptions to the empirical trends shown 
in Fig. 6 have been identified, including the widely- 
separated VLM systems 2MASS 1101-7732AB, DE- 
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Fig. 6. — fie/fj Separation (in AU) versus total system mass (in Solar masses) for known binary systems. Stellar binaries 
from Close et al. (1990); DM91; FM92; Reid and Gizis, (1997b); and Tokovinin, (1997) are shown as open squares; 
stellar-brown dwarf systems compiled by Reid et al. (2001) are shown as filled squares. The 68 binaries from Table 1 with 
measured projected separations and estimated masses are plotted as triangles; filled triangles indicate substellar primaries. 
Systems younger than 10 Myr are encircled. The dotted line indicates the maximum separation/system mass relation for 
VLM stellar and substellar binaries proposed by Burgasser et al. (2003), indicating that lower mass systems are more 
tightly bound (see also Close et al, 2003). However, three young systems (GG Tau BaBb, 2MASS 1101-7732AB and 
2MASS 1207-3932AB), and the field binary DENIS 0551-4434AB, all appear to contradict these trends. (Right) Same 
systems but this time comparing binding energy (— iSfjmd = GM1M2 / a) to total system mass. As first pointed out in Close 
et al. (2003), the widest VLM field binaries are 10-20 times more tightly bound than the widest stellar binaries, with the 
singular exception of DENIS 0551-4434AB. On the other hand, the three young VLM systems GG Tau BaBb, 2MASS 
1 101-7732AB and 2MASS 1207-3932AB are much more weakly bound. 



NIS 0551-4434AB and GG Tau BaBb. In addition, the 
extremely low mass (Mjot 0.03 Mq) brown dwarf pair 
2MASS J12073346-3932549AB (Chauvin et al., 2004, 
2005a; hereafter 2MASS 1207-3932AB), identified in the 
8 Myr TW Hydrae moving group (Gizis, 2002), falls well 
outside the mass/A„a2: limits outlined above. Such "excep- 
tions" have called into question whether current empirical 
separation limits are representative of VLM systems in gen- 
eral, and can be considered robust constraints for formation 
models; or if wide binaries are a normal (if rare) mode 
of VLM binary formation. These questions remain under 
debate. 

2.4.2 Do Evolution or Environment Play a Role in VLM 
Binary Properties? 

That three of the four weakly bound VLM systems are 
in young (<10 Myr), low density associations may be an 
important clue to their formation and existence, and encour- 
ages a closer examination of the multiplicity properties of 
such objects in general. The shaded histograms in Figs. 2 
and 3 delineate the separation and mass ratio distributions, 
respectively, of the 8 binaries in Table 1 that are members of 
clusters or associations younger than 10 Myr These distri- 
butions, although based on small number statistics, are nev- 
ertheless compelling. Young systems show a much broader 
range of separations, spanning 0.04 < A < 240 AU, with 
25tg^% (2/8) having A > 20 AU (as compared to 5tt% 



of older VLM systems). The mass ratio distribution is also 
quite flat, with a statistically significant shortfall in the rel- 
ative number of q > 0.8 binaries (25^9®% versus 8ltg%). 
Assuming that the older field sources predominately orig- 
inate from young clusters (Lada and Lada, 2003), these 
differences suggest an evolution of VLM binary properties 
over a timescale of 5-10 Myr. 

However, care must be taken when interpreting these 
data, as selection effects can distort the underlying distri- 
butions. Because the youngest brown dwarfs are still quite 
warm and luminous, imaging surveys in young clusters can 
generally probe much smaller masses — and hence smaller 
mass ratios — than equivalent surveys of older clusters or 
in the field. In addition, with the exception of some nearby 
moving groups (e.g., TW Hydra, Ursa Major), most of the 
youngest clusters lie at larger distances, so closely separated 
systems (A < 10 AU) cannot be generally resolved through 
direct imaging. This biases young samples against the close 
separations typical of field binaries. So in fact there may be 
many more closely-separated young VLM pairs, or many 
more widely- separated, small mass ratio older VLM pairs, 
than currently known. 

What about older VLM members of the Galactic thick 
disk and halo? Unfortunately, current imaging searches for 
companions to low-mass subdwarfs are not yet capable of 
detecting substellar companions directly, and radial veloc- 
ity surveys of the necessary frequency are not yet complete. 
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Gizis and Reid (2000) imaged nine VLM metal-poor (M 
subdwarf) primaries with HST, and found that none had 
companions down to the hydrogen-burning Umit. This sam- 
ple has been extended to a total of 28 M subdwarfs within 
60 parsecs, but all appear single {Riaz and Gizis, in prepara- 
tion). Taken at face value, this result {fun < 6%) suggests 
that halo VLM doubles with separations in the range 5-100 
AU are rarer than those in the disk population. However, 
given the danger of unknown selection biases, the possibil- 
ity of metallicity effects, and the still small numbers of the 
empirical sample, this result should be taken with caution. 

The current data also support the possibility that environ- 
ment may play a role in the multiple properties of VLM sys- 
tems. The three young, widely-separated binaries discussed 
above all reside in loose associations that have average stel- 
lar densities of 0.01-1 pc^-* (e.g., Luhman, 2004; Mamajek, 
2005), too low for stellar encounters to have a significant 
disruptive effect (Weinberg et al, 1987). This is in contrast 
to high-density star formation regions such as Orion, where 
average densities of lO" pc'^ (Hillenbrand, 1997) are suf- 
ficient for stellar encounters to disrupt ~10 AU VLM bi- 
naries over a ~10 Myr timescale (Weinberg et at., 1987; 
Burgasser et ah, 2003). The influence of stellar density has 
been cited for observed differences in multiphcity among 
solar-mass stars in various clusters (e.g., Ghez et al, 1993; 
Scally et al., 1999; Patience and Duchene, 2001; Lada and 
Lada, 2003; also see chapter by Duchene et al.), so differ- 
ences among VLM binaries should not be surprising. This 
scenario can also explain the paucity of wide binaries in the 
field. Dense embedded clusters, in which wide binaries can 
be easily disrupted (cf., Kroupa, 1995a,b,c) contribute per- 
haps 70-80% of the stars in the Galaxy (Lada and Lada, 
2003). The few wide systems created in less dense clus- 
ters or associations would therefore comprise a negUgible 
fraction of all VLM binaries in the field (cf., Kroupa and 
Bouvier, 2003). This scenario is compelling, but requires 
better statistics to be tested sufficiently. 

3. CONFRONTING THE MODELS 

With a full analysis of the empirical properties of VLM 
multiple systems in hand, we now examine how the pre- 
dictions of current star and brown dwarf formation theories 
compare. Detailed discussion on the current modelling ef- 
forts are provided in the chapters of Ballesteros-Paredes et 
al.. Bate et al., Goodwin et al., Klein et al. and Whitworth et 
al. Comparison of formation theories with the general prop- 
erties (mass function, disk fraction, etc.) of low mass stars 
and brown dwarfs are provided in the chapters of Duchene 
et al. and Luhman et al. Here we focus primarily on the 
predictions for VLM multiphcity. 

3.1 Fragmentation and Dynamical Evolution 

Undoubtedly, gravitational contraction of dense cores in 
molecular clouds provides the fundamental building blocks 
for stellar and substeUar objects. However, the details of 



the contraction and subsequent evolution of the cores re- 
main critical details under considerable debate. This is par- 
ticularly the case for VLM star and brown dwarfs whose 
masses are significantly below the Jeans mass (~1 M©), 
and as such cannot be formed efficiently through basic con- 
traction scenarios (e.g., Shu et al, 1987). The inclusion 
of additional physics, such as magnetic field effects (Boss, 
2001, 2002, 2004) and turbulent fragmentation, has brought 
some resolution to this problem, and has enabled a new gen- 
eration of VLM formation models. 

Turbulent fragmentation (Henriksen, 1986, 1991; Lar- 
son, 1992; Elmegreen, 1997, 1999, 2000), in which gas 
flows collide, are compressed and form gravitationally un- 
stable clumps, has pushed the fragmentation mass limit 
down to the effective opacity hmit, of order 0.01 M© (Bate 
et al, 2002). Boyd and Whitworth, (2005) have mod- 
elled the turbulent fragmentation of two dimensional sheets 
and found a protostellar mass distribution that extends to 
~ 0.003 Mq. Padoan andNordlund (2004) and Padoan et 
al. (2005) have studied three dimensional turbulent frag- 
mentation of a molecular cloud using an Adaptive Mesh 
Refinement code, and are also capable of producing cores 
as small as -^0.003 Mq. In these studies, no predictions are 
made on the overall multiphcity of the protostars. How- 
ever, fragmented cores naturally lead to the creation of 
gravitationally-bound, high-order multiple systems, as con- 
firmed in multiplicity studies of Class and I protostars 
(Haisch et al, 2002, 2004; Reipurth et al, 2002, 2004; 
Duchene et al., 2004; see chapter by Duchene et al.), and 
therefore provide a natural framework for the creation of 
VLM multiple systems. 

However, it is well known that N-body groups are gen- 
erally dynamically unstable, and dynamical scattering will 
dissolve such systems in a few crossing times (~10^ yr), 
preferentially removing the lowest-mass members (e.g., 
Kroupa et al, 1999). The scattering of low-mass bodies 
will also limit the accretion of gas and dust onto initially 
substellar cores, which would otherwise build up to stellar 
masses. These ideas have led to the so-called "ejection" 
model for brown dwarf formation (Reipurth and Clarke, 
2001), in which brown dwarfs (and presumably VLM stars) 
are simply stellar embryos ejected from their nascent cloud. 
This model has received a great deal of attention recently, 
as its qualitative multiplicity predictions - a small fraction 
of multiples and a preference for strongly bound binaries 
(close separations and near-unity mass ratios) - appear to 
fall in line with observational results. 

The most comprehensive simulations of this scenario, in- 
corporating both Smoothed Particle Hydrodynamics (SPH) 
modelling for fragmentation and accretion and N-body sim- 
ulations for dynamic interactions, have been produced by 
M. Bate and collaborators (Bate etal, 2002, 2003; Bate and 
Bonnell, 2005), and are described in detail in the chapter by 
Bate et al. Their original simulation of a 50 Mq cloud pro- 
duced only one brown dwarf-brown dwarf binary system, 
still accreting and dynamically unstable at the end of the 
simulation, implying a VLM binary fraction of <5%. It was 
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immediately recognized that this fraction may be too low 
when compared to observations (Close et al, 2003). Later 
simulations {Bate and Bonnell, 2005) found that higher 
VLM binary fractions (up to 8%) were possible in denser 
clouds. The highest density simulation also produced sta- 
ble wide (> 60 AU) VLM binary systems when low-mass 
cores were ejected in the same direction and became bound. 
It is important to note that the two wide young VLM sys- 
tems currently known are, on the contrary, associated with 
low-density associations. 

While the Bate et al. simulations have provided a great 
leap forward in the modelling of low mass star formation, 
their relevance to the observed properties of VLM bina- 
ries are hindered by necessary computational approxima- 
tions. First, sink particles encompassing all bound gas 
within 5 AU are used when densities exceed 10~^^ g cm~^. 
This approximation rules out any binaries more closely sep- 
arated than this hmit, encompassing a majority of VLM 
systems (see §2.2.2). Second, a softened Newtonian po- 
tential is employed below separations of 5 AU (down to 
1 AU), which enhances the disruption of binary pairs with 
smaller separations (Delgado-Donate et al., 2004). Again, 
as the peak of the observed VLM binary separation distri- 
bution falls within this range, it is possible that the Bate 
et al. simulations underpredict the number of VLM binary 
systems. Because the simulations are computationally ex- 
pensive, only one simulation is undertaken for a given set 
of initial conditions, resulting in poor statistics. In addi- 
tion, the simulations are allowed to run for a limited time 
(~0.3 Myr), so long term evolution of unstable multiples is 
left unresolved. 

More recent SPH + N-body simulations have attempted 
to tackle these issues by reducing the scale of the simula- 
tion. Stadiesby Delgado-Donate et al. (2004) md Goodwin 
et al. (2004a,b) have focused on smaller clouds 5 Mq) 
and have performed multiple simulations to improve statis- 
tical results. The Delgado-Donate et al. simulations were 
based on the same format as the Bate et al. work and pro- 
ceeded in two steps; first an SPH + N-body simulation of 
the gas and sink particles was conducted for -^0.5 Myr, fol- 
lowed by an N-body simulation of the resulting protostellar 
cores for a subsequent 10 Myr. This allows an examina- 
tion of both early fragmentation and accretion on the for- 
mation and disruption of bound systems, and the dynamical 
relaxation of high-order multiples over time. While brown 
dwarfs were frequently found in multiple systems contain- 
ing more massive stellar components, particularly at early 
times (~ 1 Myr), none of the simulations produced purely 
VLM binaries, again indicating a disagreement between 
theory (or at least the modelling of the theory) and observa- 
tions. A strong trend of binary fraction with primary mass 
is found, although this trend is perhaps too strong (underes- 
timating VLM multiplicity and overestimating stellar multi- 
plicity). The SPH simulations of Goodwin et al. (2004a,b), 
which tested variations of the cloud's initial turbulent en- 
ergy spectrum, also failed to produce any VLM binaries 
within 0.3 Myr. In retrospect, both sets of simulations may 



be hindered by their use of 5 AU sink particles and softened 
Newtonian potentials, and both groups have intentions to 
address these limitations (M. Bate, private conmiunication). 

Pure N-body simulations have focused on the dynami- 
cal evolution of small-N clusters of protostars, and (because 
they are less computationally intensive) have generally pro- 
duced more robust statistical predictions for VLM multi- 
ples than SPH simulations. Sterzik and Durisen (2003) 
simulated the dynamical interactions of closely-separated, 
small-N clusters and were able to broadly reproduce the 
empirical trends, including an increasing binary fraction 
and median separation with increasing primary mass (cf. 
Fig. 6), a brown dwarf binary fraction of ~10%, and a me- 
dian VLM binary separation of 3 AU. Umbreit et al. (2005) 
studied the decay of widely- separated accreting triple sys- 
tems (incorporating momentum transfer with N-body dy- 
namics) and found that VLM systems hundreds of AU apart 
were efficiently hardened to a distribution of that peaks at 
3 AU, with a long tail to wider separations. These simu- 
lations predict few very tight brown dwarf binary systems, 
although this may be because dissipative forces were not in- 
cluded. One drawback to both of these studies is that they 
do not take into account interactions with the larger star- 
forming environment, which appear to be important in the 
SPH simulations of, e.g.. Bate etal. There are plans to study 
these effects in detail (S. Umbreit, private communication). 

In short, dynamical simulations appear to reproduce 
many of the observed properties of VLM binaries, both in 
terms of quantitative results (binary fraction and separation 
distribution) and overall trends (mass dependence on binary 
fraction and mean separations). SPH + N-body simulations, 
on the other hand, generally underpredict the number of 
VLM binaries, and the lack of statistics makes the assess- 
ment of other multiple properties difficult to verify. The 
shortcomings of SPH simulations are likely related to the 
use of large sink particles. Decreasing the size of these sink 
particles, and the imposed smooth potential for close inter- 
actions, should be a priority. 

3.2 Other Formation Mechanisms 

For completeness, we briefly touch upon two other 
modes of star formation that may be relevant to the creation 
of VLM multiples. Disk fragmentation can occur when 
gravitational instabilities in massive circumstellar disks 
form, either through dynamical interactions with a pass- 
ing bare star or another disk, or spontaneously through tidal 
or spiral instabilities. Most disk fragmentation simulations 
use SPH codes to test the outcomes of different encounter 
geometries, with results depending largely on the alignment 
of the angular momentum axes of the interacting pair. The 
simulations of Lin et al. (1998); Watkins et al. (1998a,b); 
and Boss (2000) have all successfully produced substellar 
mass objects through this process (the simulations of Bate 
et al. have also produced protostellar cores through disk 
interactions). However, the disk mass necessary to produce 
such objects is nearly 0.1 Mq, and is hence unlikely around 
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a VLM primary. Therefore, while the disk fragmentation 
scenario appears quite capable of producing single brown 
dwarfs from disks around massive stars, it does Uttle to 
explain the production of VLM binary systems. 

Another VLM formation mechanism recently explored 
by Whitworth and Zinnecker (2004) is photo-evaporation. 
This process occurs when a substantial prestellar core (a 
few 0.1 M0) is compressed and stripped by the ionizing ra- 
diation front of a nearby massive O or B star. Whitworth 
and Zinnecker (2004) do not discuss binary formation ex- 
phcitly, but is possible in principle if the initial core was 
fragmented. This scenario also requires the presence of 
massive young stars, making it appropriate for high-mass 
star formation environments such as Orion, but not for low- 
mass environments such as Taurus or Cha I. Therefore, pho- 
toevaporation cannot be a universal mechanism for VLM 
multiple formation. 

4. FUTURE OBSERVATIONAL DIRECTIONS 

Despite the the large assemblage of VLM binaries now 
in place (Table 1), it should be clear that the search for 
VLM binaries should continue, particularly by broadening 
the multiphcity parameter space sampled. As such, search 
efforts should focus on low mass ratio systems {q < 0.5), 
particularly in the field; very tight systems (A < 3 AU); 
and very wide systems (A > 150 AU) with moderate to 
low mass ratios (g < 0.8). 

High resolution imaging will remain an important tool in 
the discovery and characterization of VLM binaries, partic- 
ularly with the implementation of laser guide star (LGS) AO 
systems on 5-10 m class telescopes (e.g., Palomar, Keck, 
VLT). LGS AO greatly increases the number of VLM sys- 
tems that are accessible from the ground. Ground-based 
AO enables the examination of larger samples in the field 
and in nearby moving groups and young star forming re- 
gions; and the ability to astro metrically monitor systems on 
decadal time periods, long after HST is decommissioned. 
Future studies combining AO imaging with spectroscopy 
will permit refined characterization of VLM binary compo- 
nents; note that most of the systems listed in Table 1 lack 
resolved spectroscopy. AO plus coronagraphy, the latter 
used successfully to identify several VLM companions to 
nearby, more massive stars (e.g., Oppenheimer et al., 2001, 
Lowrance et al., 2005) will facihtate the detection of low 
mass ratio systems around VLM primaries, probing well 
into the so-called "planetary mass" regime. 

For the tightest binaries, high resolution spectroscopy re- 
mains an important tool for search and characterization. Ef- 
forts thus far have been largely conducted at optical wave- 
lengths. While suitable for young brown dwarfs with M 
spectral types, optical spectroscopy becomes increasingly 
limited for L dwarfs, T dwarfs and cooler objects which are 
extremely faint at these wavelengths. Hence, future stud- 
ies should focus their efforts using high-throughput, high- 
resolution infrared spectrographs (e.g., Simon and Prato, 
2004). Short- and long-term spectroscopic monitoring cam- 



paigns of VLM samples should be pursued to identify suf- 
ficiently complete samples and to determine systemic prop- 
erties. Of the few RV variable VLM binary candidates iden- 
tified to date {Guenther and Wuchterl, 2003; Kenyon et al, 
2005; Joergens, 2006), most have only 2-4 epochs of ob- 
servation, and parameters such as separation, mass ratio, 
etc., remain largely unknown. Combining astrometric mon- 
itoring with spectroscopic monitoring for closely-separated 
resolved systems (e.g., Gliese 569Bab; Zapatero Osorio et 
al., 2004) will permit precise orbital solutions, leading to 
component mass and semimajor axis measurements, and 
enabling the examination of other multiphcity properties 
such as eccentricity distributions and spin/orbit angular mo- 
mentum alignment. 

Tight binaries can also be probed by searches for eclips- 
ing systems. For substellar objects, this is a particularly 
powerful technique, as the near-constancy of evolved (i.e., 
field) brown dwarf radii over a broad range of masses (Bur- 
rows and Liebert, 1993) implies that eclipse depths for 
edge-on geometries depend only on the relative fluxes of 
the components, while grazing transits can span a larger 
range of inclinations for a given separation. To date, only 
one eclipsing substellar system has been identified in the 
~1 Myr ONC, 2MASS J0535218-054608 (Stassun et al, 
2006). To the best of our knowledge no large surveys for 
eclipsing field VLM binaries have been undertaken. While 
eclipsing systems will likely be rare, the success and sci- 
entific yield of transiting extrasolar planet searches (e.g., 
Charbonnaeu et al, 2000) should inspire dedicated pro- 
grams in this direction. 

Interferometric observations can also probe tighter bi- 
naries than direct imaging, encouraging studies in this di- 
rection. Current facihties (e.g., Palomar, Keck, VLT) are 
limited in sensitivity, however; only the closest mid-type 
M dwarfs have been observed thus far (Lane et al, 2001a; 
Segransan et al., 2003). Increasing the throughput of these 
systems, or making use of future space-based facihties (e.g., 
SIM, TPF-1), may eventually make interferometry a viable 
observational method in the VLM regime. 

For widely- separated VLM companions, the most exten- 
sive limits to date arise from the shallow, wide-field sur- 
veys from which most of these objects were identified (e.g., 
2MASS, DENIS and SDSS). Only a few dedicated wide- 
field programs are now underway (Billeres et al, 2005; 
Allen et al. in preparation). Deep, but not necessarily 
high resolution imaging surveys around large samples of 
VLM primaries would provide better constraints on the fre- 
quency and properties of such systems. Such surveys will 
benefit from proper motion analysis and component spec- 
troscopy, allowing bona-fide systems to be extracted from 
the vast number of unrelated projected doubles. Searches 
for wide companions to young nearby stars have identified 
a few very low mass objects (e.g., Chauvin et al, 2005b; 
Neuhauser et al., 2005), and the case of 2MASS 1207- 
3934AB proves that widely separated low mass compan- 
ions can exist around VLM primaries. Future searches for 
equivalent systems, particularly in the field, will test the ve- 
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racity of the apparent wide- separation desert. 

Finally, careful selection of binary search samples 
should be of high priority. Current imaging and spectro- 
scopic field samples are largely based on compilations from 
magnitude-limited surveys, and are therefore inherently bi- 
ased. The examination of volume-limited VLM samples 
(e.g., Cruz et al, 2002) is necessary to eliminate these 
biases. Similarly, many cluster binary surveys fail to con- 
currently verify cluster membership, leading to contamina- 
tion issues (e.g., CFHT-Pl-18; Martin et al, 1998, 2000a). 
Studies have begun to address this (e.g., Luhman, 2004), 
but more work is needed. Finally, given the suggestion of 
age and/or environmental effects in binary properties, com- 
parison of large, complete samples for several clusters of 
different ages will probe the origins of multiplicity proper- 
ties and over what timescales VLM binaries evolve. 
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Table 1 

Known Very Low Mass Binaries 



Estimated Estimated Association 

Source Name Separation Spectral Types Masses q Period Age or Note Ref. 

(mas) (AU) (Mq) (Mq) (yr) (Myr) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



Cha Hq8 

2MASS J0253202+271333AB 
2MASS J0952219-192431AB 
LHS 292AB 

2MASS J21 13029-100941AB 



PPl 15AB 




0.03" 


2MASS J0535218-054608AB 




0.04*' 


2MASS J15344984-2952274AB 


65 


0.9 


GJ 569BC 


103 


0.90* 


GJ lOOlBC 


87 


1.0 


LP 349-25AB 


125 


1.3 


SDSS J092615.38+584720.9AB 


70 


1.4 


GJ417BC 


70 


1.5 


2MASS J09201 22+35 1742AB 


70 


1.5 


2MASS J2252107-173014AB 


140 


1.9 


2MASS J1847034+552243AB 


82 


1.9 


2MASS J0652307+471034AB 


170 


2.0 


DENIS PJ035726.9-441730AB 


98 


2.2 


HD 130948BC 


134 


2.4 


SDSS J042348.57-041403.5AB 


164 


2.5 


2MASS J0746425+200032AB 


220 


2.5" 


c IndiBC 


732 


2.6 


2MASS J1430436+291541AB 


88 


2.6 


2MASS J1728114+394859AB 


131 


2.7 


LP 213-68AB 


122 


2.8 


LHS 2397aAB 


207 


3.0 


LSPM 1735+2634AB 


290 


3.2 


LHS 1070BC 


446 


3.4* 


2MASS J0856479+223518AB'= 


98 


3.4 


2MASS J1017075+130839AB 


104 


3.4 


SDSS 2335583-001304AB 


57 


3.5 


2MASS J1600054+170832AB 


57 


3.5 


LP415-20AB 


119 


3.6 


2MASS J12255432-2739466AB 


282 


3.8 


SDSS 1153417.05+161546. lAB 


106 


3.8 


SDSS J102109.69-030420.1AB 


160 


3.9 


2MASS J1426316+155701AB 


152 


4.0 


2MASS J2140293+162518AB 


155 


4.0 


2MASS J15530228+I532369AB 


340 


4.4 


2MASS J 1239272+551 537 AB 


211 


4.5 


2MASS J22()6228-2047()5AB 


168 


4.5 


2MASS J0850359+105716AB 


160 


4.7 


2MASS J1750129+442404AB 


158 


4.9 


USco-109AB"= 


34 


4.9 


2MASS J2101154+175658AB 


234 


5.4 


Kelu-lAB 


291 


5.4 


2MASS J0429184-312356AB 


531 


5.8 


2MASS J0147328-495448AB 


190 


5.8 


2MASS J2I52260+093757AB 


250 


6.0 


MHO Tau SAB 


44 


6.2 


DENIS J122815.2-154733AB 


275 


6.4" 


DENIS I100428.3-114648AB 


146 


6.8 


2MASS I2147436+143131AB 


322 


7.0 


DENIS J185950.9-370632AB 


60 


7.7 


2MASS J 1 3 1 1391+803222AB 


267 


7.7 


IPMBD 29AB 


58 


7.8 


2MASS J1146345+223053AB 


290 


7.9 


CFHT-P1-12AB 


62 


8.3 


2MASS Jl 127534+741 107 AB 


246 


8.4 


2MASS J1449378+235537AB 


134 


8.5 


LP475-855AB 


294 


8.5 


DENIS I020529.0-115925AB 


510 


9.2 



M6.5 + [M6.5:] 


0.070 








M8 + [M8:] 


0.092 


0.092: 


1: 




M7 + [M7:] 


0.098 


0.098: 


1: 




M7 + [M7:] 


0.098 


0.09S: 


1: 




M6+ |M6:] 


0.085 


0.085: 


1: 




M7 + [M8:] 


0.070"" 


0.060 


0.86 


0.0159^ 


Mfi 5 + FMfi 51 


0.054^ 


0.034^ 


0.63 


0.0268*^ 




0.035 


0.035 


1.00 


4 




0.071^ 


0.054^ 


0.76 


2.4* 


T 4 S + IT 4 SI 


0.068 


0.068 


1.00 


4 


MS + TMQl 


0.090 


0.085 


0.94 


4 


T4 5 + [T4 5] 


0.050 


0.050 


1.00 


7 


T 4 S -(- n 61 


0.073 


0.070 


0.96 


7 


L6.5 + [T:] 


0.068 


0.068 


1.00 




L6 + [T2] 


0.070 


0.060 


0.86 


10 


M7 4- [M7 SI 

IVl / I [_iVl / . J 1 


0.098 


0.094 


0.96 


3 


T S + II (S 'SI 


0.075 


0.071 


0.95 


10 


MQ s + n 1 SI 


0.085 


0.080 


0.91 


11 


L4 + [L4] 


0.070 


0.060 


0.86 


14 


L7 + T2 


0.060 


0.050 


0.83 


16 


LO + LI. 5 


0.085^ 


0.066^ 


0.78 


11^ 


Tl + T6 


0.045 


0.027 


0.60 


22 


1 9 J. FT 9-1 


076 


0.075 


0.99 


15 


I 7 4- n SI 


0.069 


0.066 


0.96 


16 


MS + 11 01 


0.092 


0.084 


0.91 


15 


MS + II 7 SI 

IVIO [i-. / - J J 


0.090 


0.068 


0.76 


18 




0.082 


0.074 


0.90 


14 


MS S + FMQl 


0.070^ 


0.068* 


0.97 


16* 


L3: + [L:] 


0.071 


0.064 


0.90 


24 


L2 + [L2] 


0.076 


0.076 


1.00 


23 


LI: + [L4:] 


0.079 


0.074 


0.94 


24 


LI + [L3] 


0.078 


0.075 


0.96 


23 


M7 -4- TMQ SI 


0.095 


0.079 


0.83 


22 


T6 + [T8] 


0.033 


0.024 


0.73 


43 


T1 S + TTS SI 


0.050 


0.040 


0.80 


35 


Tl + T5 


0.060 


0.050 


0.83 


33 


MS s + n 11 


0.088 


0.076 


0.86 


27 


MS s + ri 91 


0.092 


0.078 


0.85 


27 


T7 4- rT7i 


0.040 


0.030 


0.75 


49 


L5 + [L5] 


0.071 


0.071 


1.00 


35 


MS + [MS] 


0.092 


0.091 


0.99 


31 


L6 + [L8] 


0.050 


0.040 


0.80 


39 


M7.5 + [LO] 


0.095 


0.084 


0.88 


36 


M6 + [M7.5] 


0.070 


0.040 


0.57 


46 


L7 + [L8] 


0.068 


0.065 


0.96 


49 


L2 + [L4] 


0.060 


0.055 


0.92 


52 


M7.5 + [LI] 


0.094 


0.079 


0.84 


48 


MS + [M9] 


0.0S6 


o.oso 


0.93 


47 


L6: + [L6:] 


0.069 


0.069 


1.00 


55 


M6 + [M6.5] 


0.100 


0.070 


0.70 


53 


L6 + [L6] 


0.065^ 


0.065^ 


1.00 


44^ 


LO: + [L2:] 


0.080 


0.076 


0.95 


63 


MB + [LO] 


0.084 


0.078 


0.93 


65 


LO + [L3] 


0.0S4 


0.076 


0.90 


76 


MS.5 + [M9] 


0.0S9 


0.0S7 


0.98 


72 


LI + [L4] 


0.045 


0.038 


0.84 


106 


L3 + [L4] 


0.055 


0.055 


1.00 


94 


MS + [L4] 


0.054 


0.038 


0.70 


111 


MS + [M9] 


0.092 


0.0S7 


0.95 


80 


LO + [L3] 


0.084 


0.075 


0.89 


88 


M7.5 + [M9.5] 


0.091 


0.080 


0.88 


85 


L7 + [L7] 


0.070 


0.070 


1.00 


105 



2 ChaI;RV 37 





SB2 


8; 42 




SB2 


8; 43 




SB2 


8,28; 76 




SB2 


28; 42 


120 


Pleiades; SB2 


1; 47,60 


1 


Orion; SB2, EB 


79 
5; 38 


300 


Ursa Major; triple 


2; 3,33,75,49 






25- 35 36 52 

31; 42 
69; 71 




triple 


4; 39,40 
7; 39,69,78 
32; 58,59 
23; 43 
78; 43 
4,13 




triple 


6; 26,40 

68; 43,70,71 


300 




4 7 1 7- 90 41 fit 71 


1300 


triple 


16; 40 
4; 43,86? 
4,13; 39 




triple 


17,53 

10; 36,42,82 
51' 83 




quadruple 


18; 74 
4; 43 
4; 43,86 
4' 81 

4,13; 39 


625 


Hyades 


9' 42 
5; 38,77 
15 

69; 70,72 
17; 42 
17" 42 
69; 73 
4,13; 39 
17; 42 
7; 41,52,70 
9; 42 


5 


UpSco 


29; 45,65 
4,13; 39 
24; 41,52 
23,78; 43 
78 
78 


2 


Taurus 


55; 56 

11; 41,64,71 

4 

4,13; 42 


5 


R-CrA 


20; 57 
17; 42 


120 


Pleiades 


14; 47 
7,12; 41,52 


120 


Pleiades 


14; 47,62 
17; 42 
4,13; 39 


625 


Hyades 


9; 42 




poss. triple 


12; 52,41 
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Table 1 — Continued 











E.stimated 




Estimated 




Association 




Source Name 


Separation 


Spectral Types 


Masses 


9 


Period 


Age 


or Note 


Ref. 




(mas) 


(AU) 




(Mq) 


(Mq) 




(yr) 


(Myr) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


USCO-66AB 


70 


10.2 


M6 + [M6] 


0.070 


0.070 


1.00 


120 


5 


UpSco 


29; 45,65 


2MASS J17072343-0558249AB 


950 


10.4 


M9 + L3 


0.090 


0.060 


0.67 


125 






67 


GJ 337CD 


530 


10.9 


L8 + [T:] 


0.055 


0.055 


1.00 


150 




quadruple 


30; 50,67 


2MASS J0915341+042204AB 


730 


11.0 


L7 + 1L7] 


0.070 


0.070 


1.00 


138 






78 


IPMBD 25AB 


94 


12.6 


M7 + 1L4] 


0.063 


0.039 


0.62 


200 


120 


Pleiades 


14; 47 


DENIS J144137.3-094559AB 


420 


14.3 


LI + [Li] 


0.072 


0.072 


1. 00 


200 




triple 


4,48; 39,80 


2MASS J2331016-040618AB 


573 


15.0 


M8 + [L7] 


0.093 


0.067 


0.72 


200 




triple 


4,13,17; 42,49 


USCO-55AB 


122 


17.7 


M5.5 + [M6] 


0.100 


0.070 


0.70 


250 


5 


Up See 


29; 45,65 


CFHT-P1-18AB 


330 


34.6 


M8 + M8 


0.090 


0.090 


1.00 


680 






4; 19 


DENIS J220002.0-303832.9AB 


1090 


38.2 


M8 + L0 


0.085 


0.083 


0.98 


800 






66 


2MASS J1207334-393254AB 


776 


41.1 


M8.5 + L: 


0.024 


0.004 


0.17 


2 250 


8 


TWHyd 


22,63; 34 


DENIS J055I46.0-4434I2.2AB 


2200 


220.0 


M8.5 + LO 


0.085 


0.079 


0.93 


11 500 






27 


2MASS J11011926-7732383AB 


1440 


241.9 


M7 + M8 


0.050 


0.025 


0.50 


20 000 


2 


Chal 


21; 44 



Note. — Uncertain values are indicated by colons. Additional column information: (1) name of binary; (2) angular separation in mas; (3) projected separation (A) 
in AU, or semimajor axis of orbit as noted; (4) spectral types of binary components; for sources without resolved spectroscopy, primary spectral type is for combined 
light data, secondary spectral type is estimated from photometric flux ratios (as indicated by braclcets); (5) estimated primary mass in Mq, talcen as the average of the 
reported mass ranges; masses determined from orbital dynamics are indicated; (6) estimated secondary mass in M©, taken as the average of the reported mass ranges; 
masses determined from orbital dynamics are indicated; (7) q = M2/M1, as reported or calculated from columns [5]-[6]; (8) estimated orbital period in yr, assuming 
circular orbit with semimajor axis a = 1 .26A (FM92); sources with period measurements from orbital measurements are indicated; (9) estimated age in Myr of binary if 
member of a moving group or association, or companion to a age-dated star; (10) additional notes, including cluster association; (11) references as given below; discovery 
references are listed first, foflowed by references for additional data (spectral types, distance measiu*ements/estimates, orbital measiu*ements) separated by a semicolon. 

^Parameters derived or estimated from orbital motion measurements. 

^Parameters for 2MASS J0535218-054608AB based on both spectroscopic orbit and eclipsing light curve; see Stassun et al. (2006). 
° Candidate binary. 

References. — (1) Basri and Martin (1999); (2) Martin et al. (2000b); (3) Lane et al. (2001b); (4) Bouy et al. (2003); (5) Burgasser et al. (2003); (6) Potter et 
al. (2002); (7) Reid el al. (2001); (8) Reid et al. (2002b); (9) Siegler et al. (2003); (IQ) Freed et al. (2003); (11) Martin et al. (1999); (12) Koerner et al. (1999); (13) 
G/ZK et al. (2003); (14) Marli'n el al. (2003); (15) Liu et al. (in preparation); (16) McCaughrean et al. (2004); (17) Close et al. (2003); (18) Leinert et al. (2001); (19) 
Martm et al. (2000a); (20) Bouy et al. (2004b); (21) Luhman (2004); (22) Chauvin et al. (2004); (23) Siegler et al. (2005); (24) Liu and Leggett (2005); (25) Colimowski 
et al. (2004); (26) Goto et al. (2002); (27) Billeres et al. (2005); (28) Guenther and Wuchterl (2003); (29) Kraus el al. (2005); (30) Burgasser et al. (2005a); (31) Forveille 
et al. (2005); (32) Reid et al. (2006); (33) Kmwortky et al. (2001); (34) Mamajek (2005); (35) Leggett et al. (2002); (36) van Altena et al. (1995); (37) Joergens (2006); 
(38) Tinney et al. (2003); (39) Kirkpatrick et al. (2000); (40) Ferryman (1997); (41) Dahn et al. (2002); (42) Gizis et al. (2000b); (43) Cruz et al. (2003); (44) Whittet et 
al. (1997); (45) de Zeeuw et al. (1999); (46) Kenyan et al. (1994); (47) Fercival et al. (2005); (48) Stephens et al. (2001); (50) Wilson et al. (2001); (51) Law et al. (2006); 
(52) Kirkpatrick et al. (1999); (53) Gizis et al. (2000a); (54) Delfosse et al. (1997); (55) White et al. (in preparation); (56) Briceno et al. (1998); (57) Casey et al. (1998); 
(58) Kendall et al. (2004); (59) McGovern (2005); (60) Martin et al. (1996); (61) Bouy et al. (2004a); (62) Staujfer et al. (1998); (63) Chauvin et al. (2005a); (64) 
Brandner et al. , 2004; (65) Ardila et al. (2000); (66) Burgasser and McElwain (2006); (67) McElwain and Burgasser (in preparation); (68) Burgas.ser et al. (2005b), (69) 
Burga.wrel al. (in preparation); (70) Vrha et al. (2004); (71) Geballe el al. (2002); (11) McLean el al. (2003); (73) Burgasser et al. (2002); (74) Leinert el al. (2000); 
(75) Zapatero-Osorio et al. (2004); (76) Dahn et al. (1986); (77) Burgasser et al. (1999); (78) Reid et al. (in preparation); (79) Stassun et al. (2006); (80) Seifahrt et 
al. (2005); (81) H. Bouy, private communication (2005); (82) Tinney (1996); (83) Lepine and Shara (2005) 
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